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Abstract 

We treat hadron-hadron collisions where the final state is kinematically of the 
kind associated with double-pomeron-exchange (DPE) and has large trans- 
verse momentum jets. We show that in addition to the conventional factorized 
(FDPE) contribution, there is a non-factorized (NDPE) contribution which 
has no pomeron beam jet. Within a simple model we compute DPE-two-jet 
total and differential cross sections at Tevatron energy scales, and show that 
the NDPE contribution is dominant. 
When two hadrons collide at high energy, there is about a 10% chance that one of the 
incoming hadrons will survive into the final state, while losing only a small fraction of 
its energy. Such events are called diffractive. The object that is exchanged in diffractive 
processes is called the "pomeron". A primary characteristic is that cross sections with 
pomeron exchange are approximately independent of energy at high energy; the pomeron 
has spin close to unity [|^. Moreover, there has been growing experimental evidence that 
within diffractive events, hard QCD processes can also occur, i.e., processes with a large 
momentum transfer. Such "diffractive hard scattering" has been been seen not only in 
hadron-hadron scattering by the UAl and UA8 experiments PQ], but also in electron- 
proton scattering at the HI @| and ZEUS |^ experiments, both in the deep-inelastic regime 
and in photoproduction. 

An interesting class of diffractive hadron-hadron collisions, is where both hadrons survive 
unscathed, but leave a remnant system in the central region of final-state rapidity. Such 
events are called double-pomeron exchange (DPE) events In effect, one can try to 

consider such processes as pomeron-pomeron collisions. Our purpose in this paper is to 
explore the properties of jet production in DPE events, and in particular the breakdown 
of hard-scattering factorization. Jet production by DPE has been reported in the UAl 
detector. 

The simplest model for diffractive hard scattering is due to Ingelman and Schlein 0. 
They treat each exchanged pomeron as a particle, in the same manner as one obtains pho- 
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toproduction cross sections from electron-hadron scattering at small momentum transfer. 
Diffractive hard scattering is then obtained by the use of parton densities in the pomeron. 



The model is used in much phenomenology Hfl, and has been applied to DPE []T0|. Inter 



estingly, as Collins, Frankfurt and Strikman (CFS) |lll| , p!2| explained, the Ingelman-Schlein 



model is not valid in QCD, even though the model forms a useful phenomenological bench- 



mark [|13|. CFS also predict a specific signature of the breakdown of factorization, that at 
the leading twist level all the momentum of the pomeron can go into the hard scattering — 
the process is "super-hard" or "lossless". UA8 data gives experimental support to this 
result. 

Our present paper shows how the mechanism of Collins, Frankfurt and Strikman applied 
to DPE at the level of lowest-order Feynman graphs provides a striking mechanism for the 
breakdown of factorization, and we turn it into a quantitative, if crude, modelQ The model 



is in effect a version of the Low-Nussinov-Gunion-Soper model |T^, and the same model 
was used by Berera and Soper [|I^] to understand properties of the pomeron's structure 
function. We will calculate the cross section for jet production in DPE. The important 
free parameter in our model is an overall normalization which can be determined from 



elastic scattering. Pumplin has done some calculations in the same model, but without 
the strongly dominant gluon-gluon subprocess. We call the model non-factorizing double- 
pomeron-exchange (NDPE), to contrast it with the Ingelman-Schlein model applied to DPE, 
which we call factorized double-pomeron-exchange (FDPE). 

The dramatic feature of our model is that the pomerons have no beam jets; the final state 
is then exceptionally clean, because it consists of the two isolated, diffracted hadrons, the 
high transverse-momentum jets, and nothing else. Not only are such processes theoretically 
interesting in their own right, but they have advantages for certain kinds of new particle 
searches, provided the cross section is high enough, because a lot of the background event is 
no longer present. Production processes of heavy flavors and Higgs are two such examples. 
But studies of hard scattering in DPE should surely start with the simplest processes, jet 
production. 



Schafer, Nachtmann and Schopf |]18| have characterized the final states we calculate as 
"exclusive" production , with the Ingelman-Schlein process being a corresponding "inclusive" 
process. 

Although there has been a significant amount of work on hard processes with DPE, 
almost all of it has concerned production of heavy quarks and Higgs bosons — see, for 



example, Refs. |T^|T^JT^. Moreover, much of it has used the Ingelman-Schlein method, the 



FDPE process, which, as we will see, is much smaller than the NDPE process. The only 
work on the DPE-to-jets process that we know of is by Pumplin |T^. He used the same 
model as we did, but restricted his attention to the production of quark jets. As we will 
show in Sect. |VI A| , production of quark jets is very much smaller than the production of 



gluon jets, because it has a zero at zero momentum transfer. On the other hand, Pumplin 
paid more attention to good modeling of the proton's elastic form factor. 

Donnachie and Landshoff [00] have developed an extensive phenomenology of a two-gluon 
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model for the pomeron. Bialas and Landshoff |T9[ used this model to calculate Higgs pro- 
duction by DPE. This model has many similarities to ours, and it includes better modeling 
of the non-perturbative part of the process. In particular, it has a more precise treatment 
of the pomeron-proton vertex and of the pomeron trajectory. 

However, there is one important principle that our model establishes. This is that the 
exclusive processes of DPE to jets is leading twist. This is quite non-trivial, since in proving 
factorization, there are some very non-trivial cancellations |2T|]. Indeed some of our graphs 
are a power law larger than the final answer. The proof of the necessary cancellation relies 
on Ward identities applied to the hard part, as we demonstrate around Eq. (|T3p below. To 
show that general principles do not imply some other cancellation, it is important to have 
a complete, consistent and gauge-invariant model. This is provided by our model, which 
consists of all the lowest-order Feynman graphs that follow from a certain toy Lagrangian. 

Our work also shows that one only has to be concerned about this kind of cancellation in 
the hard scattering part of the graphs. It also indicates what principles need to be applied 
in a general process. Without taking the cancellations into account, the calculation of the 
Feynman graphs is much less efficient than it need be. Of course, once these results have 
been established, one is free to do more precise modeling of the bound state. 

Lu and Milana |^ have also calculated Higgs production by DPE, but they did not 
realize that there is a factorization-breaking mechanism. Hence their cross section is too 
small, by orders of magnitude. 



I. KINEMATICS 

The process we are interested in is 

A + B^A' + B' + 2 jets. (1) 

The hadrons A and B lose tiny fractions Xa and of their respective longitudinal momenta, 
and they acquire transverse momenta Qi and Q2. (This defines a diffractive regime, and in 
Regge theory would lead to an expectation of the dominance of double pomeron exchange 
(DPE) — Fig. |l].) The jets carry large momenta of magnitude Et in the plane perpendicular 
to the collision axis with azimuthal angle 0. (This defines a hard-scattering regime.) The 
small transfer of longitudinal momentum to the hard process implies large rapidity gaps 
between the jets and the two outgoing hadrons. In what follows, we will generically denote 
the hadronic scale by m. We are interested in the kinematic region where m, Qi, and Q2 
are of a typical hadronic scale (less than about 1 GeV), while Et is much greater than this 
scale, but much less than ^/s. 

To understand the asymptotics of Feynman graphs for our process, it is convenient to 
use light-cone components (+, — ; _L). The components of momenta of the hadrons in Fig. |l| 
are 
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FIG. 1. Double pomeron exchange (DPE) to two jets. 



2 (1 — X )v2s 

We use bold-face type to indicate two-dimensional transverse vectors. For the jets, we ignore 
terms of relative order ^ 1 to get: 

Pi = [aXaJ -,bxb^ -•,Et cos (f),ET sin (f) 



P2 = (bxa^, axb^^■, -Et cos 0, -Et sin 0^ , (3) 



where it is convenient to define 



with 



1 + 

2 ' 

6 = 1 - a, (4) 



K = —^. (5) 

X(xX})S 



We are interested in the limit where Xa^Xb 0, to give a Reggc- style limit, and where 
Et oo, to give a hard scattering limit. We hold k, Qi, and Q2 fixed, to correspond to a 
fixed angle for the hard scattering and fixed transverse momenta for the outgoing hadrons. 
In this limit, s necessarily goes to infinity, since it is proportional to E^/xaXb- 
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\^ ^ Pi' 





FIG. 2. Our model of the non-factorizing Double-Pomeron-Exchange (NDPE) amplitude with 
two gluon jets produced. 

II. MODEL 



Our method is that of Berera and Soper ||T6| . We model the process by the lowest order 
Feynman graphs that are appropriate. The pomerons are replaced by two-gluon exchanges, 
while instead of true bound states, we model the hadrons by elementary color-singlet scalar 
particles that we will call "mesons" and that are coupled to scalar quarks by a 0^ coupling. 

We normalize the coupling of our mesons so that we reproduce the measured value for 
cross section for small angle elastic scattering of protons, when we use two-gluon exchange 
for elastic scattering. This is essentially the Low-Nussinov-Gunion-Soper model |T^, but 
with some simplifications that should not change the relative sizes of the cross sections too 
much. (We will address the adequacy of our approximations in Sect. [VI1| .) 

Our model for jet production by double pomeron exchange (DPE) is therefore given by 
Fig. 0. Two gluons couple to each of the diffracted hadrons: this is the minimum number of 
gluons that can couple to a color-singlet. The top and bottom bubbles imply a sum over all 
one loop graphs with a quark loop — Fig. ^. One pair of gluons scatters to make the two 
jets, as in Fig. ^. Note that since the mesons are color-singlet, the jet pair is color singlet 
and hence Fig. ^ does not include the graph with an s-channel gluon. Also in Fig. |^ the 
graph with a two-gluon two-quark contact interaction is omitted, since, as we will see after 
Eq. @, this graph gives a non-leading contribution in the Regge limit. . 

For our purposes, the most essential feature of the pomeron is that it gives cross sections 
that are approximately independent of the center of mass energy ^/s, when ^/s is made 
large with all transverse momenta fixed in size. Such behavior is given by exchanges of 
spin-1 fields — see Thus, in a perturbative model, the lines coupling bubbles with very 
different rapidities must all be gluons, exactly as in Fig. |[ Exchange of a quark-antiquark 
pair instead of a gluon pair would give an amplitude suppressed by a power of s. Moreover, 
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FIG. 3. Model for two-gluon form factor of meson. There are three other graphs identical to 
the above, except with the arrows on the quark lines reversed. 




FIG. 4. Graphs for g + g ^ jet + jet, when the jets are gluon jets, given that the jet pair is 
color-singlet. Similar graphs apply for production of a quark and antiquark jet. 

for the same reason, the two-gluon form factors in Fig. ^ must be between the meson pairs 
{A, A') and {B,B'), rather than between the pairs {A,B') and {B,A'). 



III. CALCULATION OF AMPLITUDE 

To compute the amplitude in Fig. H, one first performs contour integrals for the lon- 
gitudinal components k~ and of gluon loop momentum, while making the appropriate 
approximations to get the leading power behavior in the limit we are considering. In the 



leading region one finds ||23| 



Pi, ^ {^JCaVs, m'^/y/s, rnj , (6) 



P2' ~ {^"^/y/s, Xb\/s, m 
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The following approximations can now be made: 



• Neglect k~ in all but the top bubble. 

• Neglect in all but the bottom bubble. 

• Take the dominant component (+ times — ) in the product of the currents at the end 
of each gluon line. 

In addition, the integral over all but one light-cone component of momentum can be easily 
done for the loop integral inside each form factor. The integral over the transverse momen- 
tum k still links the different parts of the integral. Since individual graphs for the hard 
scattering Fig. ^ are a power E^/m^ bigger than the sum, we cannot yet neglect the loop 
transverse momentum in the hard scattering. 

The approximations can be represented as follows. Let the exact formula for Fig. ^ be 

J [271^ k^ pf, P2, 

where pi> = k + p^—pA' and = —k + ps—PB', while F represents the 2-gluon form factor 
of the meson. Fig. ^, and A represents the hard scattering subgraph. Fig. ^. The symbol 
/ labels the final state of the hard scattering. Then the approximations we have made are 
equivalent to replacing this by 

/ ^F++{iO,k-,0),pA,PA')F-ii-k+,0,0),PB,PB')Ai-+-J;pr,p2')^^^. 
J [Air) — K P2' 

(8) 

(The bold-face type indicates two-dimensional transverse vectors.) We now can see that, 
as stated earlier, the graph where there is a two-gluon-two-quark vertex in Fig. |^ gives a 
zero contribution in Eq. (H). After doing the contour integrals, we find that the momentum 
components of the gluons entering the hard scattering are 



Within the hard scattering amplitude A, we will be able to replace these vectors by 
(xqa/s, 0, 0) and (0, Xbi/i, 0), but only after taking care of the cancellations that reduce 
it by a factor of order m^/E^. The gluons pi' and p2' will be referred to as the incoming 
gluons for the hard scattering. 

The expressions we obtain for scattering amplitudes can conveniently be written in terms 
of the following quantity obtained from the meson form factor^ 



'Compare the results for the impact factor in the theory of elastic scattering in [23|. 
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iV2-lJ(k)V(Q + k)2^o 



da {Gf{a, aQ) - Gf{a, k + aQ)} 



(10) 



where / labels the incoming hadrons {A ot B), a and b are color indices, A^'c = 3 is the 
number of colors, and where 



Gf (a,v) = -g'G''6' 



2^2 c-ab 



a(l — a) 



(2vr) 



In 



4A(a) 



4A(a) 



4A(a) 



(11) 



Here, 



A(a) = - a(l - a)M^ 



(12) 



and now m represents the quark mass. To obtain the integral in Eq. (|10|) for I = Awe start 
with the sum of graphs in Fig. ^ (plus the three graphs with the quark lines reversed). In 
accordance with the approximations that give Eq. (^), we take the + component for each of 
the vector couplings, we set the = 0, and we apply 



dk- 



(13) 



Finally we divide by P^. 

This gives the integral in Eq. ( p!OD for the case I = A, and a similar procedure is apphed 
for c/B- The prefactors in Eq. (|l^) are adjusted so as to absorb the color factors and the 
denominators of the gluon propagators in a convenient way. 

There is now one trick needed to simplify the calculation of the amplitude Fig. |^, to 
exhibit the cancellations we mentioned earlier. We have used A{ij,, f) to represent the 
amplitude for the hard scattering part. Fig. ^ Here, fi and u are the Lorentz indices for the 
incoming gluons pi' and p2', and / is a label for the jet part of the final state. As explained 
above, we need only the component A{—,+;f). However, individual Feynman graphs for A 
are a factor of order / bigger than the final answer. Following standard methods, we 
add terms that are zero by the QCD Ward identities and obtain: 



Ai-,+-J)=Ai-,+-J) 



Pi'^,A{fx,+; f) A{-,u-J)p2'u , Pv^,A{iJ,,iy; f)p2^ 



+ 



Pi' P2' 

After dropping terms that are non-leading by a power, we obtain 



PT'P2' 



(14) 



Ai-,+-J) 



PviAjiJ; f)p2\ 
PvPt 



(15) 



where the latin indices i and j refer to the transverse components of vectors. In Eq. ([15|) , the 
factor pviP2'il {Pi'P2') is of order m'^/E^, which demonstrates the previously mentioned can- 
cellation. Within the hard scattering amplitude A{i,j] /), it is now correct to approximate 
pi> by {xa\fsj2, 0, 0) and p2' by (0, Xb\[sj2, 0). 
For the amplitude we now find that 
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M = -^-^ f -^^^(k,-Qi)^B(k,Q2)6,(k-Qi)e,(-k-Q2)^(^,j;/) 



XaXb J {2tt) 

-iBij{Qi,Q2)A{i,j;f) 



Xa^b 

(16) 



Here, the "polarization" vectors are defined as 



e.(k) = (17) 



It will be convenient to define 



M = XaXbM. (18) 



This scaled amplitude has the property of being independent of Xa, Xb, and Et in the 
combined Regge and hard scattering limit that we are considering. 

We note for later use that when the transverse momenta of the two outgoing hadrons 
are zero, the hadronic tensor Bij in Eq. (|16D is invariant under rotations about the collision 
axis, so that we can write 

Bij{0,0) = 2SijP^{0). (19) 

(We will see later that f3 is related to the meson-pomeron coupling in our model; its definition 
here agrees with Eq. ( p8D below.) It will also be convenient to extract the couplings and to 
define 

B, - (20) 

In Eq. (^) A{i,j] /) is the amplitude for an incoming gluon pair in linear polarization 
state i,j to go to the outgoing final state /. The incoming gluon pair for the hard scattering 
will be in an overall color singlet state since the remaining portion of the final state, the two 
mesons, is color singlet. We will refer to A{i,j]f) and -Bij(Qi,Q2) as the hard scattering 
amplitude and the hadronic amplitude, respectively. Both amplitudes are 2x2 tensors in 
the transverse space of gluon polarizations. 



IV. CROSS SECTION 

We now give formulae for the cross section. In terms of the light-cone momentum frac- 
tions Xa and Xb, 

da \M\^K^ 

^"^^ d^Qid^ Cl2dEld(t)jdxadxb ~ 2i77r8E4yr^ ^ ' 

with implicit sums over final state color, flavor and spin. Next let us change to the rapidity 
variables of the jet partons in the final state. These are defined as 
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so that 



_ 1 



where y„ = yi — ?/2- Note that depends only on In terms of this variable and 

y+ = {yi +y2)/2 we have 

da m-'n' ^24) 



d^Qid'^Q2dE^d(f)jdy-dy+ 

The squared amplitude can be expressed in terms of the hadronic and hard amplitudes 
as follows: 

\M\^ = B*^H,^kiBki, (25) 
where Bij is defined in Eq. ([l^) and 

H,,kM=Y.A*{i,j-J)A{k,l-f). (26) 

/ 

In the hard amplitude A{i,j; /), i and j are gluon polarization indices while / generically 
represents any final parton pair state. The sum Z)/ is over all spin, flavor and color states 
for final-state partons of given momenta. 

From Eqs. (p|)-(|25|) at fixed k, we see that the differential cross section is leading twist 
and has the appropriate large s behavior to approximate an amplitude with pomeron ex- 
change. That is, when Ex gets large and Xa and Xh get small, the amplitude Ai is constant, 
and the cross section goes like 1/E^. 

The easiest way to see that this is a leading twist contribution is to verify that the 
same power laws are obtained from photon exchange. Such a process is obtained when one 
replaces the hadrons in reaction (|l]) by charged particles, and the exchanged pomerons in 
Fig. |l| by single photons. Normally, when one has a short-distance subprocess initiated by 
elementary particles like photons, and one replaces the photons by a composite objects like 
pomerons, one loses a power of the scale E^ of the hard scattering. But this result fails when 
a diffractive requirement is imposed on the final state, as explained by Collins, Frankfurt 
and Strikman . Our model is perhaps the simplest example of a failure of hard-scattering 
factorization. We have exhibited all the graphs at a certain order of perturbation theory for 
a particular final-state, and there is no possibility of cancellation. The factorization theorem 
is an order- by-order result; to get factorization one has a cancellation between the graphs 
that we calculate and some other graphs with different final states. 



V. NORMALIZATION FROM ELASTIC SCATTERING 

The parameters of our model are: the mass m of the scalar quark, the mass M of the 
mesons, the quark-meson coupling G, and the gauge coupling g. We suppose that the most 
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important unknown is the normalization of our model for the non-perturbative physics. So 
we set all the masses to a typical hadronic scale: = = 0.1 GeV^. As for the couplings, 
we have a factor in the hard scattering amplitude A and a factor g'^G'^ in the hadronic 
part. The g"^ in the hard amplitude should be given by the usual running coupling at a scale 
of order Et, and we now determine a suitable numerical value for g^C^ from applying our 
model to elastic scattering. This model is effectively the Low-Nussinov model [p!5| . 
We write the elastic-scattering amplitude for our model in the Regge-like form 

M,i{t) = sP^t). (27) 

Here, the pomeron-hadron coupling is given by a transverse integral, 

(3'{t) = -\g'j |^^A(k,Q)^B(k,Q), (28) 

where Q is the momentum transfer, so that t = — Q^, and g{k, Q) is defined in Eq. (|10|). 
We also define 

- (29) 

The elastic scattering cross section is then 

=^ = ,ScsM^. (30) 

We can now determine gG by equating the left-hand side of this equation to the pp-elastic 
cross section at some particular t = and at some appropriate value of s. From Eq. (|30D 
we can reexpress the DPE amplitude |A^p (defined in Eqs. (|T6|) - (^)) as follows: 

|A7p = 16.f^') ^llElE^^ (31) 

where the unknown couplings have dropped out. 

This equation, (|3TD, is particularly simple for forward scattering of A and 5, when we 
set all the momentum transfers to zero: Qi = Q2 = io = ^- Then we use Eq. (|19D to give 



|A<(0, 0)|2 = 647r 5,,5uiH,,ku (32) 

which relates the DPE-to-jets amplitude to an elastic cross section and a hard scattering 
amplitude The elastic cross section is observable, and hard scattering amplitudes are per- 
turbatively calculable. 



VI. NUMERICAL RESULTS 

In this section we compute numerical values for the DPE-to-jets cross section. 
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A. Zero in forward quark-jet production 



Now, there are two hard subprocesses: gg —>■ gg and gg — > qq. We first show that 
the gg qq subprocess gives a zero in its contribution to the DPE-to-jets process when 
the transverse momentum transfers from the hadrons are zero, a result previously obtained 
by Pumplin Since the gg gg gives no such zero, it must be by far the dominant 

subprocess, and Pumplin's estimates for the cross section should be much lower than the 
true values. 

By a straightforward perturbative calculation, one can verify that for the quark ampli- 
tude, gg qq, 



where Aj is the helicity of quark i. This property implies, from Eq. (|32D , that the quark 
amplitude vanishes when the hadrons A and B are forward scattered. That this same 
cancellation does not occur for the gluon amplitude can be verified from Eq. (|A5|) in the 
appendix. 

The cancellation in the case of quark jet production can be understood by examining the 
helicity amplitudes given by Gastmans and Wu |]2^ for e~^e~ — >• 77. (Notice that because 
our final state is color- singlet, only the abelian parts of the graphs are relevant.) 

Hence our cross section is dominated by the production of gluon jets. All the polarization 
combinations for the gg — > gg amplitude are given in the appendix, in Eq. ( |A5|) . Special 
attention is given in deriving these amplitudes with incoming linearly polarized gluons. 
In terms of the amplitudes in the appendix, the final state sum in Eq. ( pB] ) requires a 
polarization sum and a factor of 8 for color: 



For comparison purposes, we have also computed (a) the inclusive two-jet cross section 
(Fig. ^) (i.e., without a diffractive requirement: A + B — > jet + jet + X), and (b) the 
result of applying the Ingelman-Schlein model to DPE, which gives a result for the process 
A + B ^ A' + B'+iet + jet -|- X. This last process we call the factorized double-pomeron- 
exchange (FDPE) two-jet cross section, and it is given by graphs like Fig. ||. 

For inclusive two-jet production, the differential cross section is. 




(33) 




(34) 



B. Comparison cross sections 



dEj^dy^dy. 



da 



Y.^'^f<^/A{Xa)Xbfb/B{Xb) 



a,b 



(35) 



where 




(36) 
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FIG. 5. Amplitude for inclusive two-jet production. 




FIG. 6. Factorized-Double-Pomeron-Excliange (FDPE) amplitude with two jets produced. 
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and -ffinci is the standard hard inclusive matrix element squared for incoming partons a, b. 
It is summed over all final parton pairs. Note that it is averaged over initial polarization 
states of the incoming partons a and b. 

For FDPE, the natural extension of the Ingelman-Schlein model gives 



da 



dEj,dy^dy^ 



2567cE^ Vl67r 



J dxp^ dxpg Bf^ {xp^ ) Bf^ {xpg ) 



Xafa/PA i^a/ XpJXbfb/Pg {Xb/ Xpg)H- 



ab 
incl 



(37) 



where is a normalization factor |T3| that we have set equal to 2/7r to give the conventions 
of Donnachie and Landshoff . We have used the following integral over the momentum 
transfer: 



Bl{xp)= r dt \(3i{t)\\] 

J —oo 



(JoXp 



2ao 



(38) 



where ap(t) = + C('t is the pomeron trajectory function. We will use yao = 4.6 mb2, 
bo = 3.8 GeV"^ ao = 1.08 and a' = 0.25 GeV-^ following m. 



C. Results 

In Fig. 1^, the differential cross section, da/dE^dy+dy_, is plotted in mb/GeV^ at y/s = 
1800 GeV, y^ = and Et = 5 GeV. We give the differential cross section for the three 
cases: inclusive jet production, for factorized DPE, and for our model of non-factorizing 
DPE. 

In the calculation of inclusive cross sections, the parton distributions in the proton are 
those of CTEQl [^. In the calculation of the FDPE cross section, the parton distributions 



in the Pomeron are given by using the CTEQ package to evolve the distributions numerically 
from initial distributions at = 4 GeV^. As initial distributions, we use an ansatz like 
that of Ingelman and Schlein . For the gluons we choose 

/,/p(x)=3(l-x), (39) 

and for the light quarks we choose 

3 

fu/pix) = fd/p{x) = -(1 - x); (40) 

we set all other parton species to zero. These initial distributions satisfy the momentum 
sum rule, and split the pomeron's momentum equally between gluons and quarks. For our 
estimates of cross sections, these parton densities are sufficiently close to fits obtained at 
HERA 1^ . In all our calculations, the distribution functions were evolved with three quark 
fiavors and A = 0.281 GeV. In the hard amplitude, the running coupling constant was 
evaluated at Et- 
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For NDPE, we normalize the couplings in our model of the hadrons from the forward 
elastic scattering cross section, as explained in Sect. 0. We use the value {da{0)/dt)ei = 
201 mb/GeV^, which is obtained by the optical theorem from a total cross section a^^t = 
62.7 mb p6|p8| , |29[| . Since the true elastic cross section depends on energy, but our model's 
cross section does not, we chose to use the measured cross section at -^i = 540 GeV. so 
that the pomeron is being sampled under approximately the same conditions in both the 
DPE-to-jets and the elastic scattering process. The NDPE cross section depends linearly on 
the elastic cross section — see Eq. (|3TD — so that one can trivially rescale our predicted cross 
sections when a different value for {da{0)/dt)ei is preferred. 

The total cross sections integrated over y_ and for Et > 5.0 GeV are, at y/s = 
1800 GeV, 

ainci(1800,5) = 2.4 mb, 
c^fdpe(1800,5) = 0.0022 mb, 
(Xndpe(1800,5) = 0.68 mb. 



(41) 



and at = 630 GeV, 



o'inci(630, 5) = 0.31 mb, 
c^fdpe(630,5) = 0.000062 mb, 
cxndpe(630,5) = 0.18 mb. 



(42) 



The integration range for the rapidities ?/+ and was restricted so that both incoming 
partons into the hard process had momentum fractions less than 0.05. This corresponds to 
a typical selection cut on the diffractive hadron for identifying pomeron exchange events. 
The same cut on the incoming partons was made for the estimate of the inclusive jet cross 
sections. If one allows the complete range of momentum fractions from to 1 for this case 
one obtains crinci(1800, 5) = 5.0 mb and crinci(630, 5) = 1.3 mb. For quark jets in the NDPE 
case, we have verified that their cross sections are at least a few orders of magnitude smaller. 
Also for the NDPE case, by a simple hand calculation of the total cross section using Eqs. 
(P^) and (^) one can confirm the order of magnitudes quoted above. From this exercise, the 
factor 4 increase at = 1800 GeV compared to i/i = 630 GeV is seen to arise primarily 
from the same factor increase in the accessible region of jet rapidity. 

We will comment on the implications of these calculations in the next section. 



VII. CONCLUSION 

We have estimated jet production in DPE processes, by a mechanism in which the whole 
of the momentum of the pomerons goes into the jets. The process has a quite dramatic 
signature: the final state consists of the two diffracted hadrons, two high-i^^ jets, and 
nothing else. It is a leading-twist contribution, not suppressed by a power of Et, despite the 
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FIG. 7. Differential cross section da / dEj^dy+dy- at v^=1800 GeV, Et=5 GeV for: top solid 
two-jet inclusive, bottom solid 2 X (NDPE), and dashed 200 X (FDPE). 
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fact that the pomeron is a composite object. This is permitted because the factorization 
theorem does not apply to diffractive processes. 

Moreover, the cross section is an order of magnitude larger or more than the obvious 
Ingelman-Schlein type process. The lack of pomeron beam jets is presumably the reason 
why we are able to get such large cross sections. 

Our model is, of course, quite primitive. However, by being a complete and consistent 
calculation in lowest order perturbation theory, it does establish an important principle: that 
one cannot treat the exchanged pomeron as a particle with associated parton densities. The 
approximations we have made are concerned only with extracting the leading power behavior 
in the appropriate kinematic limit. Since gauge-invariance causes some quite non-trivial 
cancellations before the final result is obtained, the exact gauge invariance and consistency 
of our model is critical to establishing the general principles. 

We have also a new relation, Eq. (|3^) , which relates the DPE-to-jets cross section at zero 
momentum transfer to the elastic scattering cross section. This relation is independent of 
all details of the hadronic form factor, and depends only on the use of a two-gluon-exchange 
model. 

There are at least two obvious sources of large corrections to our model: 

• Absorptive corrections due to extra exchanges of pomerons and gluons between those 
objects in our model that have very different rapidities. Recent work by Gotsman, 



Levin and Maor [30] gives one way of tackling this problem. 



Sudakov corrections that suppress the process of two gluons of low virtuality going 
into a hard process. Standard QCD methods can presumably be applied here. 



From Ref. PU|, for example, we expect that the true DPE-to-jets cross section will be an 



order of magnitude or more smaller than the values we have calculated. But remember that 
absorptive corrections will reduce all DPE cross sections, including FDPE. More precise 
modeling of the hadron form factor is also needed. But this is presumably a less important 
matter. 

We plan to return to these issues in a later paper. 

In addition, a comparison with existing data from the UAl detector would be useful. 
Since more experimental and phenomenological work is necessary to obtain correctly normal- 
ized cross sections, we will merely note here some orders of magnitude. At y/s = 630 GeV, 
Ref. PI quotes a cross section for detected DPE events of 0.3 fih, whereas Streng's M 



estimate for the same quantity is quoted as 5 /ib. This discrepancy is attributed to the 
stringent cuts imposed on the data. 

Moreover, the experiment finds that 49% of their events have one or more jets with Et 
above 5 GeV. Normally, in a hadron-hadron collision, such a jet fraction would be impossibly 
high. For example, our FDPE cross section, in Eq. (^), after reduction by absorption and 
experimental cuts, is much smaller than the data. But our NDPE cross section is tens of /xb, 
so until we apply absorption and Sudakov corrections this partial cross section is in danger 
of being larger than the full DPE cross section. 

What we may reasonably conjecture is that hard scattering cross sections in DPE may 
be surprisingly high. Although our cross sections are about an order of down from the fully 
inclusive jet cross section (with the same cuts), it should be remembered that the DPE 
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cross section (without the jet requirement) is itself much smaller than the total cross section 

iiH. 



We suggest that a qualitative understanding of the large cross sections can be obtained 
by observing that, in our model, Fig. 0, we really do not have pomerons exchanged between 
the hadrons and the jets, but only single gluons. So that we get a color singlets coupling 
to the hadrons, it is sufficient to add one extra gluon exchanged between the two diffracted 
hadrons. 

If further examination supports the sizes of cross sections we predict, then the mechanism 
we are using could be very important for all kinds of studies. For example, one can produce 



the Higgs boson [O. Although the cross section would be a lot lower than the total Higgs 



cross section, the lack of a background event could make up for the lower rate in terms of 



the usefulness of the signature, at least for certain ranges of parameters. Compare Ref. pT 
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APPENDIX: FORMULAS FOR AMPLITUDES 

We give the amplitudes for 

g{pv, ev{i)) + g{p2', e2'(j)) ^ 9{Pi, ei(«)) + g{p2, e2(/5)), (Al) 

We present them with the color factor 6cd for the final-state gluons removed. That is, we 
write the amplitude A{i,j; /) in Eq. as 

A{t,j-J) = 6,,A{z,j;a,(3). (A2) 



Note also that to get A normalized as in Eq. ([Iq) we contracted the incoming gluons' color 
indices with 5""^ /{N"^ — 1) = 5ab/^- Our notation for the polarizations (ei'(«), ei(a), etc.) 
indicates that we are using a different set of basis states for each of the four gluons. A 
gluon's polarization vector must be transverse with respect to its momentum: • p = 0. 
For the incoming gluons, we define linear polarization vectors as 

ei,(l) = e2'(l) = (0,0;l,0) 

ei,(2) = e2'(2) = (0,0;0,l), (A3) 

and for the outgoing gluons, we define polarization vectors in (/) and out (O) of the x-y 
plane as 

= 62(1) = (O,O;sin0, -COS0), 
ei{0) = 

^2(O) = W£^(l,^;0). (A4) 
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Then straightforward Feynman graph calculations give: 



^(l,l,/,/) = i(2,2,0,0) = z/^ 
i(l,l,0,0) = A(2,2,/,/)=z/^ 
i(l,2,7,0) = i(2,l,0,7) = -i/ 

A(l,2,0,/) = A(2,l,/,0) = ^/^ 



— 1 ) sin^ (j) + cos^ ( 



— 1 ) cos^ 6 + sin^ I 



- cos^ (p+T sin^ (p 
a b 

1 2.1.2 

- COS H — sm I 
b a 



A{1, 1, /, O) = A{1, 1,0,1) = -A{2, 2, /, O) = -A{2, 2, 0, 1), 

' sin (p cos 0, 



• 23 
'^4 



A(l, 2, 1, 1) = A(2, 1, /, /) = 2, 0, 0) = -i(2, 1, 0, 0) 



sm cos ( 
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